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Abstract 



We show how to "concatenate" variational principles over different 
bases into one over a single base, thereby providing a unified Lagrangian 
^ treatment of interacting systems. As an example we study a Klein- 

3 Gordon field interacting with a mesically charged particle. We employ 

our method to give a novel group-theoretic derivation of the kinetic 
stress-energy-momentum tensor density corresponding to the particle. 
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1 Introduction and Setup 

Let us recall the geometric setting of a classical variational principle (GIMMsy I 
[1998]): We are given a fibration Y ^ X, with dimX = n + 1, and we wish 
to extremize an action of the form 

SW = [ C{j'^) 



J X 

where %p : X -^ Y is a. section and C : J^Y -^ A"+^X is a specified Lagrangian 
density.^ 

One commonly encounters several (say K) such variational principles si- 
multaneously, for instance when one studies the Newtonian dynamics of a 
swarm of charged particles (in a background electromagnetic field), or the in- 
teraction between Dirac and Yang-Mills fields. In the cases cited, the relevant 
fibrations have the form Yi -^ X for the i^^ variational principle; the key point 
being that each fibration has the same base X. To combine these variational 
principles into a single principle is a straightforward matter: one builds the 
fiber product Yi x^ • • • XxYr -^ ^, and then on the first jet of this bundle one 

takes as the Lagrangian density an expression of the form £i H h Ck + -^int 

for some interaction terms £int. 

It is less clear how to deal with variational principles with disparate bases, 
that is, fibrations Yi —* Xi in which the Xi are all different. A simple example 
is a nucleon moving in a dynamic Klein-Gordon field. (Here the configuration 
bundle for the nucleon is X x M ^ M, where X is 4-dimensional spacetime and 
M is the material world line of the nucleon. The fibration for the Klein-Gordon 
field is M X X — > X, sections of which are scalar fields on spacetime.) Even 
if the bases are identical, it may be desirable to distinguish them. This is the 
case, for instance, in relativistic multiparticle systems (cf. Anderson [1967]), 
when one wants to parametrize each particle's trajectory by its own proper 
time, as opposed to a single "universal" time. 

In this context of disparate bases, one standard way to proceed is as follows. 
Construct an action functional using sections ipi : Xi ^f Yi for the z*^ bundle 
by setting 

5(^i,...,V'i^) = V / A(jV.)+ I £i„t(jVi,---,jVx). (1.1) 

Then varying these fields ipi, one obtains the Euler-Lagrange equations for the 
problem. See equation (2.1) for a specific example. 

However, while producing the Euler-Lagrange equations, this approach has 
the unsatisfactory feature of not yielding a field theory in the usual sense, in 



^ For simplicity we consider only first order theories. We also ignore technical issues and 
proceed formally. 
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which the fields are sections of a single bundle and which has a well-defined 
Lagrangian density. This or some other formalism is needed if one wishes to 
tap into the machinery of multisymplectic geometry, multimomentum maps, 
stress-energy-momentum ("SEM") tensors, and constraint theory, etc. 

To concatenate variational principles with disparate bases in such a way as 
to recapture a genuine field theory, we proceed as follows. To begin, construct 
the product bundle YiX- ■ -x Y^ -^ Xi x ■ ■ ■ x X^, which we denote F — > X for 
short. In agreement with experience we restrict attention to product sections 
of this bundle of the form tp = [tpi, . . . jipx), where each ipi is a section of 
Yi -^ Xi. With ip = {ipi, . . . , ipx) such a section, 

J V(a;) = (jVi(a;i), • • • ,fipK{xK)) 

where x = (xi, . . . , xk)- Denote by J^Y the subbundle of J^Y consisting of 
all such jets; equivalently, J^Y = J^Yi x ■ ■ ■ x J^Yk- 

Given Lagrangian densities Ci on the jet bundles J^Yi, it is simple enough 
to lift them to maps, still denoted by Ci, on the concatenated jet bundle J^Y 
by composing with projections: 

i^il){x) i-^ CiQ^ipiixi)). 

But how do we concatenate these Ci into a single Lagrangian density? Even 
ignoring interaction terms, we cannot just add the Ci as they take values in 
different spaces, viz. A"'+^Xj and so need not be forms of equal rank. The 
trick is to "suspend" the d : .PY -^ A^^^+^X^ to maps PF -^ A^+^X, where 
N = ni + ■ ■ ■ + uk, by inserting suitable tensor densities in the Ci to "even 
out" their ranks in the target. 

First, we pull Ci back via the projection X -^ X^ to an (n^ + l)-form on X. 
Second, for each i choose scalar densities 1)i of weight 1 on Xi x ■ ■ ■ Xj ■ • ■ x Xk- 
Now in Ci = Lid^^^^Xi the coefficient Li transforms as a scalar density of weight 
1 on Xj, so the coefficient in 

Ci ■= Lid''^+\ A S),rf"^+ixi A ■ • • A (f^i A ■ ■ ■ A rf"^+V 

= ±LiT)id^+^x 

will also transform as a scalar density of weight 1 on X under the subgroup 

Diff(Xi) X ■ ■ ■ X Diff(XK) C Diff(X) 

(which is sufficient for our purposes). The densities SDj are to be chosen by 
hand, depending on the precise structure of the system; see the examples in 
§§2 and 3. Thus modified, we may assemble Ci + ■ ■ ■ + Ck into a map 
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Interaction terms, which are typically defined over several of the bases Xi 
(again, see the following examples) are treated similarly. Finally, it is straight- 
forward to deal with composite situations in which some of the bases are iden- 
tical and others are not. 

Ultimately, the specific choice of the 2)j will not matter as long as 



X J X 



Lj-iJijj^ d oc — / i-j-i d X 



for each i, that is, the concatenated action reduces to the original action. 
Specifically, this means that 

/ £(j V) = SW 
Jx 

where the right hand side is given by (1.1). In particular, the Euler-Lagrange 
equations remain unaltered when the Lagrangian L is used in place of the action 
functional (1.1). 

Once we have a total Lagrangian density in hand (albeit possibly a distri- 
butional one), we may proceed in the usual fashion. Thus we may compute 
the equations of motion and various geometric objects, such as SEM tensors. 
To extract physical information from these objects, however, it will normally 
be necessary to "project" them from X to some Xi or products thereof; this 
projection is accomplished by integration over the remaining Xj. 

Rather than continuing to try to describe the procedure in generality, it is 
more instructive to illustrate it via a simple example. (It really is easier done 
than said!) 

In §2 we apply this method to a system consisting of a Klein-Gordon 
field interacting with a mesically charged particle. (Think of a pion field 
interacting with a nucleon.) Beyond illustrating concatenation, this example 
has interesting features which are worth elucidating. In particular, we study 
the SEM tensor density of this system. Its computation, following Gotay and 
Marsden [1992], is interesting in that it naturally produces the Minkowski, 
or kinetic, SEM tensor for a moving particle as a matter of course. To our 
knowledge, this SEM tensor has never been derived via a Lagrangian from first 
principles; it has always been inserted into the formalism in an ad hoc manner. 
An important point therefore is that our method is not merely a 'tidy' means 
of packaging variational principles; it is capable of providing, in an entirely 
straightforward fashion, quantities which otherwise cannot be obtained except 
in makeshift ways. 

Finally in §3 we briefiy indicate some other contexts in which our results 
should be useful. 
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2 Motion of a Mesically Charged Particle in a 
Klein— Gordon Field 

Let X be an oriented spacetime with metric G. We consider a real Klein- 
Gordon field : X — * M of mass M interacting with a particle of mass m and 
mesic charge e. The particle's trajectory in spacetime (or "placement field") 
is 2 : M — > X. The base for the system is thus X x M, the second factor being 
thought of as a time axis,^ and the configuration bundle Y is then 

(M X X) X (X X M) ^ X X M 

with coordinates (0, X") on the fiber and (x'^, A) on the base. We set z"" = 
X'^oz. 

Our presentation is based upon the excellent exposition in Chapter 8 of 
Anderson [1967], to which we refer the reader for further information. The 
action (1.1) for the system in this case is usually written 



5(0, z)= f^l (G"^'^(a;)0,^(x)0,,(x) - M2</)(x)2) v^^ 
e(f){x)\\z{X)\\6\x-z{X))d^xdX 



XxR 

- f m\\z{X)\\ dX, (2.1) 

Jr 

where the dot denotes differentiation with respect to A and ||i|| = \/—GabZ°'z''. 
Observe that the bases for the free Klein-Gordon term and the free particle 
term are different, and that the interaction term in the middle lives on the 
product of these. 

Before proceeding, there are two technical issues that need to be resolved, 
stemming from the presence of the two factors of X in the configuration bundle. 
First, note that in the leading term of S, G is regarded as living on the X 
in the base, while in the last term it evidently resides on the X in the fiber. 
It is necessary to know precisely where G lives, as this has an effect on the 
subsequent analysis: if on the base, then G is treated as a field, while if on 
the fiber it is simply thought of as a geometric object. We reconcile these two 
interpretations by taking G to be anchored to the base, and then pulling it back 
to the fiber by means of the following construction.^'^ Introduce yet another 



^ Not necessarily proper time. 

^ This is a variant of the Kuchaf method of parametrizing a classical field theory; see 
Gotay and Marsden [2008] and Castrillon Lopez, Gotay, and Marsden [2008] for details. 

** At the end of this section we will briefly examine what happens if instead we anchor 
G to the fiber. 
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factor of X in the fiber along with diffeomorphisms t] : X ^ X, viewed as 
sections oi X x X ^ X, with corresponding configuration and multivelocity 
variables 77° = X"or/ and 77"^ = d{X°' ori)/dx^, respectively. (We can, and do, 
regard the two copies of X in the fiber as identical.) We use these auxiliary 
nondynamic fields, the covariance fields, to {i) identify the copies of X in the 
fiber with that in the base, and {ii) endow the new copy of X in the fiber with 
the metric g = rj^^G with components 

9ab = G^yK af^ ft, 

where (n^a) = {''1°' n)~^- AH this is summarized in the figure below. 



Fibers, 



X X X 



Base, X : 




Klein-Gordon field, <i> 



covariance field, r] 




metric, G 



The general set up for the introduction of covariance fields. 



Second, the delta function (5^(a; — -2(A)) must be modified, as it compares 
elements x in the base with elements z{\) in the fiber. As just indicated we 
can use the covariance fields to remedy this problem as well: we need only 
write 5^[x — ?7^^(2;(A))) instead. It is sometimes convenient to replace 



5\x-r\z{\))) = 5\7]{x) - z{\)){Aei7],) 



(2.2) 



using the properties of delta functions (cf. the Appendix), where rj^ is the 



Jacobian of rj. From this we see that 5 



X — Tj 



j(A))) (z) is a scalar density on 



X (again, see the Appendix), and {ii) depends upon the spacetime derivatives 
of 7], even though this is not obvious at first glance. The reason we do not 
insist on a fixed identification of the base X with the fiber X, and instead 
allow a variable identification by means of the covariance fields, will become 
apparent in due course. 
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Remark. Analogous fields t], called covariance fields as well, are introduced 
in Gotay and Marsden [2008] and Castrillon Lopez, Gotay, and Marsden [2008] , 
but there they have a different purpose, namely, to make a field theory on a 
given background generally covariant and in doing so, they are introduced as 
dynamic fields. ♦ 

In addition to the covariance fields rj, we introduce a (positive-definite) 
metric A' on M as a nondynamic field. We suppose that K is chosen so that 
M has metric volume 1. 

With these fixes we may now concatenate the three action terms over the 
composite base X x M as S{(j),z) = j^^^Ld'^xdX, with the Lagrangian 

L(a;^ A, </), <^,^, /, i^ 77^ r/%, G,., i^) 

- [m + e(t))\\z\\5^{ri - z){deiri^). (2.3) 

Notice that the interaction term itself needs no essential modification, as the 
corresponding term in (2.1) is already an integral over X x M, but informs our 
choice of scalar density in the free particle term, viz. 5^{x — ?7~^(2;(A))), when 
we suspend the latter to X x M. We have also written this delta function in 
the form (2.2) to make it clear that L is defined pointwise. 



Remark. The choice of scalar density 2) = yK in the Klein-Gordon term 
is hardly unique; all we require is that ^^1) d\ = 1. For instance, we could 
instead take 6{X) for 2) with no essential difference. ♦ 

As evident from (2.3), the modified configuration bundle is taken to be 

Y' = Y XxX^Xx Lor(X) x^ Riem(M), 

where we abbreviate the bundle X x X -^ X by X^, Lor(X) is the bundle 
whose sections are Lorentz metrics on X and, similarly, Riem(]R) is the bundle 
whose sections are Riemannian metrics on M. However, in our approach (p and 
z are variational, while r], G and K are nondynamic fields. As per the above, 
we now regard 

Remark. Occasionally, as in Landau and Lifshitz [1979], one encounters 
what one might call "noncovariant concatenations." In the current example, 
this amounts to writing the terms in the action as integrals over X alone and 
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is effectively accomplished by imposing the coordinate condition x° = A. As 
this procedure is not covariant, it can lead to problems (Leclerc [2006]). ♦ 

We compute the Euler-Lagrange equations. Varying with respect to and 
employing (2.2), we obtain 



-MV(x) a/^GR /?r(A) - E\\ziX)\\ 6\x - r/-i(^(A))) 

- d^ (g^>,,v^=g) (x) ^/K{X) = 0. 



Integrating with respect to A, using the fact that voli^( 
this reduces to the Klein-Gordon equation 



VV^cf) + M'cj) = -p 
where V denotes the G-covariant derivative and 

p(x) = e(-G)-5 f ||i(A)|| S^{x - ri-\z{\))) d\ 



1 , and rearranging, 
(2.4) 



R 



is the source density. 

Similarly, varying with respect to z and employing (2.2) yield 



1. 
dz"" L 



{m + e(P{x))\\zi\)\\5\x-ri-\zi\))) 
d 



+ 



dX 



(rn + e^ix))'-^!^^^^^;^ 



0. 



Carrying out the differentiation and then integrating over X, some manipula- 
tions give 



d 
d\ 



{m + e<j>{r^-\z{X)))) 



g,,{z{X))z\X) 

mm 



-6K' 



'aM^''i<m)uim 



+ {m + e(f){r]-\z{X)))) 



,Mmz\X)z%X) 

2l|i(A)|| 



(2.5) 



To give insight into these equations, note that in the special case when {X, G) 
is Minkowski spacetime, rj = Idx, and A is taken to be proper time along the 
particle's world line, these equations simplify in a global Lorentz frame to 



d 
dX 



(m + e(f){z{X)))za{X) = -£0,a(2;(A)). 
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This is the mesic analogue of the Lorentz force law in electrodynamics. 

Neither K, the G^jy, nor the 7]°" have field equations, since they are not 
variational. Thus one is free to assign them whatever values one wishes in 
(2.4) and (2.5). Often, however, one has specific values of G and rj in mind, 
e.g., the given spacetime metric for G and Idx for rj. 

Turning now to the SEM tensor, let Diffc(^) x Diffc(IR) (that is, the group 
of diffeomorphisms that are the identity outside a compact set) act on the 
modified configuration bundle Y' according to 

(ax/)- (x. A, 0, z] T], G, K) = (o-(x), /(A), 0, z; r], a^G, f^K) . 

(We assume that all diffeomorphisms are positively oriented.) The Lagrangian 
density C. = L d*x dX is then visibly equivariant with respect to the induced 
action on J^Y', that is,^ 

m^ X /) • j\<p, z- 77, G, ir) = (a X f)X{3\4>. z; V, G, K)) . 

We may thus use equation (3.12) in Gotay and Marsden [1992] to compute 
the b- dimensional SEM tensor density 



'T^l n-A 

rr \ ■'■ V ■'■ V 

' ' A -^ 4 



of the interacting system (where x^ = A).^ Integrating over A and raising an 
index, we project out the spacetime SEM tensor density: 



Jr 
= i^"" + (m + e0)e'^^ (2.6) 



where 



t^^ = -- [pG^'^G"^ - G^'^G"^)0,„0,^ + G^'^'M^cjAV^ 
is the canonical SEM tensor density of the (free) Klein-Gordon field and 

e^'^(x) = k\k\ 1 4^)^ 5'{x - r\zW)) dx 



^ Even though the pointwise action of DifF(X) on the fiber of the "covariance bundle" 
X X X —> X is taken to be trivial, its action on sections thereof is not: a -rj = rjo a~^ . Thus 
the identification of the factor of X in the base with that in the fiber can fluctuate, which 
is why we allow rj to be variable in the first place. 

^ Using the product metric G Q) K on X x M., one could also compute T via the Hilbert 
formula (4.2) in Gotay and Marsden [1992]. See also Leclerc [2006]. 
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is the Minkowski tensor density. {mQ^^ is then the kinetic SEM tensor den- 
sity). As well, we compute T^;^ = = T'^4. Finally, we find that when 
integrated over X, T^4 is effectively the Klein-Gordon action: 

Remark. The kinetic SEM tensor density is a familiar object in microscopic 
continuum mechanics, cf. Chapter 8 of Anderson [1967] and §33 of Landau 
and Lifshitz [1979]. Minkowski [1908] originally introduced it in flat-spacetime 
electrodynamics in order to recover the continuity equation Ti^'^ ^u = in view 
of the fact that T^m,!^ t^ when currents are present. In the continuum limit of 
a noninteracting clutch of particles, Q^" goes over to the SEM tensor density 
for a perfect fluid as in §§9.1-2 of Anderson [1967]. It is interesting that in this 
limit, the inflnite time integrals in the kinetic SEM tensor density disappear 
and one is left with a local tensor density. 

To our knowledge, ours is the flrst genuine derivation of the Minkowski 
tensor density from flrst principles in a variational context, once again illus- 
trating the power of multisymplectic geometry in classical fleld theory and in 
particular, the usefulness of having a concatenated theory for which one can 
make use of concepts such as the SEM tensor. 

As we have deflned it, the Minkowski tensor density depends upon the 
covariance flelds as well as the particle placement fleld. However, note that 
when Tj = Idx, © reduces to the more familiar expression 

Remark. Suppose we focus solely on the particle dynamics so that the (orig- 
inal) conflguration bundle is X x M — > M. The corresponding Lagrangian 
density — ?7i||i(A)||(iA is Diffc(]R)-covariant, and so we may compute the cor- 
responding SEM scalar density as in Example a. Interlude II of GIMmsy II 
[2004]. We obtain T = —E, the energy of the particle, which vanishes as 
the Lagrangian is time reparametrization-invariant. (This is reflected by the 
vanishing of T^4 in the 5-dimensional context when there is no Klein-Gordon 
fleld.) Thus only when the spacetime X is part of the base of the variational 
principle do we encounter the kinetic SEM tensor density; it does not appear 
in standard particle dynamics per se. 

To reiterate, even in the absence of other flelds, our technique yields yet 
another (5-dimensional!) treatment of the relativistic free particle that has 
the advantage of automatically incorporating the Minkowski tensor. ♦ 
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Remark. Note that the term e(j)Q^^ in (2.6) arises from the interaction 
of with the mesically charged particle. This term has no analogue in the 
electrodynamics of particles; there we get simply 

Charged strings behave similarly, as we show in §3A (cf. equation (3.1)). This 
can be traced to the fact that the electromagnetic field is a covector, while the 
Klein-Gordon field is a scalar. ♦ 

The SEM tensor density T is symmetric. It is also divergence- free, as can be 
seen from general principles (cf. Proposition 5 in Gotay and Marsden [1992]). 
One may verify this directly, via a long calculation. 

We end with a discussion of an alternate treatment of this system. 

Suppose we consider the physical metric as a geometric object g on the 
fiber as opposed to a field on spacetime. Then we would define G = rj*g with 
components G^y = rj"-^rf'vgah- Proceeding as in the above, the Lagrangian 
density would be 

L(x^A,0,0,^,/,i^r/^r/%,ir) 

= \ («:^«:^^°V,M0,. - MV') v^ (det r/,) ^/K 
— {m + e(t))\\z\\5^{ri — z){detri^) 

where a/— G = ^/^(det?]*) and ||i|| = ^J—gab{z)z"■z^. 

Gomputing the SEM tensor density in this formulation, we obtain T^'j, = 
and the other components as before. That the spacetime components vanish is 
actually a consequence of the generalized Hilbert formula (3.13) in Gotay and 
Marsden [1992], since the nondynamic fields r] and K do not transform under 
Diff(X). (In the original formulation, the nondynamic metric G on X does 
transform under the spacetime diffeomorphism group with the result that (2.6) 
is nonzero.) The difference between this SEM tensor density and the previous 
one stems from: {i) the spacetime metric no longer being regarded as a field, 
so that it cannot contribute to the energy, momentum, and stress content of 
the system, and {ii) the subtly different manners in which the covariance fields 
appear in the two formulations. 

That one can encounter several SEM tensor densities for the 'same' system 
may seem surprising, but is unavoidable and can also be regarded as different 
"packaging" of the same information. What the SEM tensor density turns 
out to be depends upon what the fields are, whether they are dynamic, and 
precisely how they appear in the Lagrangian. And even the size of the SEM 
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tensor density depends upon how the system is formulated! For instance, for 
something as simple as a relativistic free particle, we can have a 1 x 1 SEM 
tensor density (which vanishes identically) — as noted in a previous remark, or 
a 5 X 5 SEM tensor density (which does or doesn't, depending on where the 
spacetime metric is anchored). And in the latter case, the 5x5 object reduces 
to the 4x4 Minkowski tensor density! Thus how the system is formulated 
plays a substantial role insofar as how various quantities, and in particular the 
SEM tensor density, are to be understood. 

3 Further Examples and Outlook 

To conclude we briefly mention some other systems for which our techniques 
should prove helpful. We begin by upping the dimension of the matter from 
1 to 2, that is, we replace the particle by a string. For variety, we also replace 
the mesic interaction by an electromagnetic one. 

Charged Strings. We closely follow the exposition in §2. Let (X, G) and 

(W = 'Rx B,H = —HkQHb) be 4- and 2-dimensional Lorentzian spacetimes, 
respectively. We consider a string worldsheet in X, this being a map z : W ^ 
X. We use coordinates (x'^, A"^ = (t, a)) as coordinates on X x W. The 
configuration bundle Y is correspondingly 

A^X Xw{X xW)xx Lor(X) Xw^ot{W) ^ X x 1^. 

Assume that the string carries a charge density p: S — >■ M and interacts 
with a dynamic electromagnetic field described by a potential 1-form A on X. 
We also take the metric if on VF to be dynamic; thus we adopt the Polyakov 
approach as in Green, Schwarz, and Witten [1987]. The action for this system 
is 

S{A,z,H) = -\ [ F^''{x)F^,{x)V-G{x)d^x 

+ / A^{x)p{a)——{T,(7) 6'^{x — z{T,a)) y'HBicr) d\dT da 
Jxxw '^'T 

-^ f H^''{\)G,M\))z^{X),Az''{\),B^/^H{\)d^\ 
^ Jw 

where T is the tension. 

As in the case of the meson, we see that z takes values in X, which is the 
base for the electromagnetic field. So we need to introduce covariance fields as 
before. As well, we take the spacetime metric to reside on the factor of X in 
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the base. Finally, let K he a nondynamic Riemannian metric on W with total 
volume 1. The modified configuration bundle for the concatenated variational 
principle is then 

Y xx X^ Xw Riem(Vr) ^XxW 

and the Lagrangian reads 

+ A^ k\p -q:^S^{v- z) (det 7],) ^/Hb 

- ^if^^ G^, k\k\z\a z\b S\v - ^)(det 7],)V^H. 

Now build the product metric G(BK on X x W. Using the Hilbert formula 
and then integrating as before, we compute the 6-dimensional SEM tensor 
density as follows: the spacetime components are 

%^''' = %'^^ + Te>"' (3.1) 

where 

is the free electromagnetic SEM tensor density, and 

e'^'^ = K^^K^b / H^^z\Az\BS\x-r]~\z{a)))V^Hd^a 
Jw 

is the analogue of the Minkowski tensor density for strings. The 'extra' T^^^ 
and T^^ components are zero and, after integrating over X, the T^^ subblock 
reduces to 

Continua. Another intriguing example that we intend to pursue in future 
works is a charged elastic body, fiuid, or plasma, in which one concatenates 
a continuum with electromagnetism on a given background metric spacetime. 
Such theories will likely have significant differences with the particle and string 
examples presented above. Amongst these differences, we expect that, unlike 
mesically or electrically charged particles, continua should have well-defined 
initial value problems (see also the discussion of this point in Anderson [1967]). 
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Evidence for this can be found in works such as Kunzinger, Rein, Steinbauer, 
and Teschl [2005] and Rein [1990]. 

One other interesting aspect of a charged elastic body is the following. If 
B is the body manifold, then its motion in spacetime is determined by a map 
z : M. X B —>■ X. The main difference from our previous two examples is that 
rather than the delta functions 6^{ri{x) — z{X)) we must now use characteristic 
functions x{v~^{^(M ^ -^)))- We expect that examples such as this will be key 
players in the future development of the point of view given in this paper. 

Appendix 

Let M be a manifold with coordinates x = {x^, . . . , x"^). Here we prove that 
the delta function 6"^{x — Xq) transforms as a scalar density of weight 1. 
Let ?7 : M — > M be a diffeomorphism and / G C°°(X). On the one hand. 



(/o^)M = ifor]){x)5^ix-xo)d"'x. 
Jx 

On the other hand, by the change of variables theorem with y = ri{x), 

fivixo)) = [ f{y)S"'iy-vixo))d^y 



X 



firi{x))S^ir]ix)-ri{xo))\Jix)\d^x. 



X 



where J is the Jacobian determinant of rj. Since / is arbitrary the desired 
result follows upon comparing these two formulae. 
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